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This  study  focuses  on  gaining  fundamental  insights  into  the  sparsely  explored  area  of  reacting  multiphase 
flows  in  a  compact  heat  exchanger.  Chevron  plate  heat  exchangers  have  been  demonstrated  to  possess 
superior  thermal  performance,  scalability,  and  mixing  capability  compared  to  more  traditional  shell- 
in-tube  heat  exchangers  or  stirred  tank  batch  reactors.  This  study  explores  the  hydrodynamic  behavior 
of  gas-evolving  reacting  flows  in  chevron  plate  heat  exchangers.  Experimental  characterization  of  a  plate 
heat  exchanger/chemical  reactor  (in  multi-phase  flow  configuration)  utilizing  high-speed  video  and  axial 
pressure  measurements  was  conducted.  Existing  correlations  that  were  developed  using  air-water  flow 
in  PHEs  predicted  with  acceptable  accuracy  the  total  pressure  drop  in  the  HEX  reactor. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  concept  of  using  a  compact  heat  exchanger  as  a  continuous 
flow  chemical  reactor,  termed  as  heat  exchanger  (HEX)  reactors,  is 
well  documented  in  the  chemical  engineering  literature  as  dis¬ 
cussed  by  Anxionnaz  et  al.  [1].  The  authors  cite  high  volumetric 
heat  transfer  capacities  ranging  from  1400  to  4000  kW/m3.  This 
represents  2-3  orders-of-magnitude  gain  in  the  value  of  surface- 
area-to-volume-ratio  when  compared  to  that  of  batch  reactors. 
Among  the  many  classes  of  compact  heat  exchangers  examined 
by  these  authors,  the  gasketed  plate  heat  exchangers  (PHEs)  offer 
both  high  thermal  performance  and  ease  of  maintenance.  PHEs 
can  be  easily  disassembled  for  inspection  and  maintenance  (in  con¬ 
trast,  shell-and-tube  heat  exchangers  require  very  elaborate  main¬ 
tenance  procedures  due  to  the  complexity  of  assembly / 
disassembly  1]).  Also,  the  performance  rating  of  a  gasketed  PHE 
can  be  scaled-up  simply  by  assembling  additional  plates.  The 
advantages  of  vigorous  mixing,  high  surface  area  and  low  volume 
make  the  PHE-based  HEX  reactor  an  ideal  platform  for  continuous 
delivery/removal  of  heat  into/from  a  reacting  multi-phase  flow  [2]. 

Although  the  continuous  flow  HEX  reactor  represents  a  promis¬ 
ing  approach  to  enhancing  heat  transfer  in  applications  involving 
high  energy  density  gas-evolving  reactions,  the  hydrodynamics  of 
combined  multiphase  flow  and  chemical  reaction  within  a  PHE 
are  not  well  understood.  Therefore,  the  motivation  of  this  study 


*  Corresponding  author  at:  A236,  Building  18, 1950  5th  St.,  Wright-Patterson  AFB, 
OH  45433,  United  States.  Tel.:  +1  (425)  358  1295. 
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is  to  gain  insights  into  the  nature  of  reacting  multi-phase  flows 
in  a  PHE  serving  as  a  HEX  reactor.  The  objective  of  this  paper  is 
to  establish  an  experimental  basis  for  the  hydrodynamics  of  gas- 
evolving  reacting  flows  using  a  model  chemical  reaction,  and  to  as¬ 
sess  the  accuracy  and  suitability  of  existing  correlations  for  void 
fraction  and  pressure  drop  calculations  in  a  chevron  plate  heat  ex¬ 
changer.  The  hydrodynamic  correlations  evaluated  in  this  paper 
are  used  in  Part  II,  which  focuses  on  developing  numerical  models. 


2.  Experimental  methods  and  materials 

2A.  Chevron  plate  heat  exchangers 

Chevron  PHEs  consist  of  a  stack  of  thin  heat  transfer  plates  with 
corrugated  surfaces  in  the  shape  of  a  chevron  (Fig.  1).  These  plates 
are  separated  by  compliant  gaskets,  which,  when  the  entire  stack  is 
mechanically  compressed,  seal  the  space  in  between  a  pair  of  mat¬ 
ing  plates  to  form  channels.  The  gasket  also  serves  to  direct  the 
cold/hot  streams  to  their  respective  ports,  as  can  be  seen  in 
Fig.  1.  The  chevron  patterns  of  the  upper  and  lower  plates  that 
form  an  individual  channel  face  opposite  directions  to  form  a  tor¬ 
tuous  network  of  crisscrossing  passageways.  Fluid  streams  direc¬ 
ted  along  the  upper  and  lower  corrugations  interact  to  produce 
unique  secondary  flow  patterns  that  give  rise  to  enhanced  convec¬ 
tive  heat  transfer  [3-5].  Additional  details  of  the  experimental  set 
up  and  measurements  are  provided  in  an  earlier  study  [2]. 


0017-9310/$  -  see  front  matter  ©  2013  Elsevier  Ltd.  All  rights  reserved. 
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Nomenclature 

A 

surface  area 

X 

mass  quality 

C 

molar  concentration 

z 

axial  coordinate 

D 

diameter 

L 

length 

Greek 

M 

molecular  weight 

A 

corrugation  pitch 

P 

pressure 

two-phase  multiplier 

pP 

partial  pressure 

Ot 

void  fraction 

Re 

Reynolds  number 

chevron  angle 

T 

temperature 

p 

dimensionless  concentration 

U 

superficial  velocity 

P 

density 

V 

volume  flow  rate 

W 

plate  width  (between  gasket) 

Subscripts 

X 

Lockhart-Martinelli  parameter 

A 

acceleration 

F 

friction 

G 

gravity 

Lowercases 

SP 

single-phase 

d 

V2  Corrugation  depth 

TP 

two-phase 

f 

Darcy  friction  factor 

c 

corrugated  section 

g 

gravitational  acceleration  constant 

e 

equivalent 

k 

rate  constant 

g 

gas  phase 

k * 

dimensionless  rate  constant 

h 

hydraulic 

m" 

gas-phase  mass  flux 

l 

liquid  phase 

Ktal 

total  mass  flux 

PP 

port-to-port 

t 

time 

U 

velocity 

2.2.  Experimental  setup 

The  hydrodynamic  effects  of  a  solid-gas  reacting  flow  was  char¬ 
acterized  using  the  well-known  acid/base  neutralization  reaction 
between  acetic  acid  and  sodium  bicarbonate: 

NaHC03  +  CH3COOH  ->  CH3COONa  +  H20  +  C02  (1) 

The  reaction  in  Eq.  (1 )  is  non-toxic,  and  occurs  at  room  temperature 
and  atmospheric  pressure;  these  attributes  make  acetic  acid  and  so¬ 
dium  bicarbonate  an  ideal  model  reaction  to  study  the  hydrody¬ 
namics  of  gas-evolving  multi-phase  reacting  flows. 

The  rate  of  reaction  (and  gas  evolution)  was  controlled  by 
varying  the  concentrations  of  the  reactant  species  being  fed 
through  the  reactor,  and  the  reactor  was  run  under  adiabatic  con¬ 
ditions.  The  overall  reaction  in  Eq.  (1)  is  first-order  with  respect 


to  acetic  acid  [6]  and  is  diffusion  limited  in  aqueous  solution 
[7].  Because  of  the  high  degree  of  mixing  that  occurs  in  PHEs, 
the  concentration  of  the  reactants  is  controlled  so  that  the  major¬ 
ity  of  the  reaction  is  completed  inside  the  channel,  and  the  vol¬ 
ume  of  gas  produced  can  be  determined  directly  from  the 
molar  concentrations  of  acetic  acid  and  sodium  bicarbonate  intro¬ 
duced  into  the  reactor. 

The  HEX  reactor  is  an  Alfa  Laval  MF3  gasketed  chevron  plate 
heat  exchanger  having  the  following  dimensions,  as  detailed  in 

Table  1. 

To  aid  in  flow  visualization,  the  top  plate  of  the  channel  was  re¬ 
placed  with  a  transparent  window  fabricated  from  an  acrylic  cast¬ 
ing  resin  (Electron  Microscopy  Sciences,  Hatfield,  PA).  The  resin 
was  molded  over  one  of  the  PHE  heat  transfer  plates  to  form  the 
chevron  pattern,  so  as  to  preserve  the  channel  geometry  while 


Fig.  1.  Alpha  Laval  MF3  heat  transfer  plate  used  in  this  study  (left)  and  profile  of  corrugation  (right). 
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Table  1 

Geometric  data  of  Alfa  Laval  MF3  PHE. 


A  (mm) 

11 

De  (mm) 

4.68 

P (deg) 

30 

Lpp  (mm) 

323.6 

Lc  (mm) 

275 

W  (mm) 

94.3 

allowing  for  observation  and  video  acquisition  of  the  flow.  The 
window  thickness  is  12.7  mm. 

The  PHE  was  mounted  vertically  and  the  flow  directed 
upwards,  as  indicated  by  the  red  arrows  in  Fig.  2.  Dilute  solutions 
of  sodium  bicarbonate  and  acetic  acid  were  fed  into  the  reactor 
at  the  inlet  port,  subsequently  mixing  and  initiating  the  reaction. 
Fig.  3  shows  the  flow  loop  layout  for  the  reactor,  which  consists 
of  two  circuits  for  the  sodium  bicarbonate  solutions  that  can  be 
operated  independently  to  adjust  the  molar  ratios.  Each  of  the 


two  flow  loops  features  a  turbine  flow  meter  and  a  centrifugal 
feed  pump  (Tuthill,  Alsip,  IL)  controlled  by  a  10-turn  potentiom¬ 
eter.  The  window  is  tapped  with  5  instrumentation  ports  to 
which  diaphragm  absolute  pressure  transducers  with  a  50-psi 
ceiling  were  attached.  The  uncertainties  provided  by  the  manu¬ 
facturer  for  the  above-stated  instruments  are  tabulated  in 
Table  2. 

The  pressure  transducers  and  flow  meters  interface  with  a  dig¬ 
ital  data  acquisition  system  comprised  of  National  Instruments’ 
SCXI  series  data  acquisition  cards  (1102  for  the  pressure  transduc¬ 
ers  and  1126  for  the  flow  meters),  which  is  in  turn  controlled  by  a 
LabVIEW®  virtual  instrument  (VI)  (National  Instruments,  Austin, 
TX)  for  data  logging  and  real-time  monitoring.  High-speed  video 
was  acquired  with  a  Phantom  V4.2®  camera  (Vision  Research, 
Wayne,  NJ)  fitted  with  a  100  mm  Macrovision  lens  (Carl  Zeiss,  Ger¬ 
many).  The  camera  was  positioned  approximately  0.2  m  from  the 
front  of  the  acrylic  window. 

2.3.  Procedure 

The  effects  of  flow  rate  and  molar  concentrations  of  reacting 
species  on  the  flow  pattern  and  pressure  profile  within  the  PHE 
channel  were  evaluated.  The  chemical  solutions  were  mixed  prior 
to  conducting  the  experiments:  99%  sodium  bicarbonate  (Sigma 
Aldrich,  St.  Louis,  MO)  was  dissolved  into  a  4  L  mixture  of  distilled 
water  and  99%  propylene  glycol  (Sigma  Aldrich,  St.  Louis,  MO)  to 
obtain  a  mass  concentration  of  20g/L.  The  volume  fractions  of 
water  and  propylene  glycol  were  further  adjusted  to  vary  the  mix¬ 
ture  viscosity  and  thence  the  range  of  Reynolds  number  attainable 
by  the  pumps.  Acetic  acid  solutions  were  mixed  in  1  L  batches  from 
glacial  acetic  acid  (Sigma  Aldrich,  St.  Louis,  MO)  and  distilled  water 
at  volume  ratios  (acid: water)  between  5%  and  10%.  The  relative 
uncertainties  calculated  using  the  Kline-McClintock  method  [8] 
were  ±0.44%,  0.60%,  and  0.13%  for  the  propylene  glycol  volume 
fraction,  acid  solution  volume  fraction,  and  sodium  bicarbonate 
concentration,  respectively. 


Y  Waste  Line  Valve 


Nall  CO,  Reservoir 


Acid  Injection  Line 


#-t£> 


Acid  Bypass  Valve 


Fig.  3.  Flow  diagram  of  HEX  reactor;  P0-P4  indicate  pressure  transducer  positions. 
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Table  2 

Instrument  uncertainty  provided  by  manufacturer. 
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Instrument  name 

Manufacturer 

Model 

Uncertainty  ( ±  ) 

NaHC03  Flow  meter 

Omega  engineering  (Stamford,  CT) 

FTB2002 

3.0%  of  reading 

Acid  flow  meter 

Omega  engineering  (Stamford,  CT) 

FTB2002 

3.0%  of  reading 

Pressure  transducers 

Omega  engineering  (Stamford,  CT) 

PX409-050A5  V 

0.08%  of  reading 

Prior  to  conducting  the  experiments,  the  window  surface  was 
polished  and  coated  with  a  thin  layer  of  carmuaba  wax  to  protect 
the  acrylic  resin  from  the  acetic  acid.  The  lines,  valves,  and  flow 
meters  were  flushed  with  distilled  water  before  and  after  testing 
to  remove  contaminants.  The  following  procedure  was  adopted: 
the  bypass  valves  on  both  flow  loops  were  opened,  and  their 
respective  pumps  adjusted  until  the  desired  flow  rate  was  obtained 
for  each.  During  this  time,  the  illumination  apparatus  and  high¬ 
speed  camera  position  were  adjusted.  The  pumps  were  temporar¬ 
ily  stopped,  the  bypass  valves  closed,  and  the  waste  valve  opened; 
the  camera  and  pumps  were  then  triggered  simultaneously,  during 
which  time  the  VI  recorded  pressure  and  flow  rate  data  while  the 
camera  captured  images  at  1000  frames  per  second.  After  10  s, 
the  acid  pump  was  shut  off  and  the  sodium  bicarbonate  pump 
was  allowed  to  continue  running  for  an  additional  5  s  to  purge 
any  remaining  acid  from  the  chamber.  The  camera  was  then  ad¬ 
justed  to  a  new  axial  position,  and  the  test  repeated  twice  at  the 
same  acid  and  sodium  bicarbonate  flow  rates. 

3.  Results  and  discussion 

The  experiments  were  conducted  by  varying  the  fluid  flow  rate 
and  the  molar  ratio  between  the  reacting  species.  The  flow  rates 
ranged  from  0.031  to  0.074  L/s;  the  acidrbicarb  molar  ratio  ranged 
from  5:2  to  7:10.  The  conditions  for  each  experiment  are  tabulated 

in  appendix  A. 

3.1.  Pressure  measurements 

To  establish  that  the  duration  of  the  individual  experimental 
runs  were  sufficient  to  permit  a  statistically  steady  [9]  condition 
to  be  reached,  the  absolute  pressure  recorded  for  a  representative 
test  set  are  plotted  in  Fig.  4  as  a  function  of  time  elapsed  from  the 
test  start.  The  pressure  and  flow  rates  undergo  a  brief  startup  tran¬ 
sient  once  the  acid  pump  is  started,  which  quickly  plateaus  after  an 
elapsed  time  of  5  s,  and  remains  steady  until  the  acid  pump  is 
stopped.  This  behavior  was  observed  for  all  molar  ratios  and  flow 
rates  considered,  so  time  averaged  data  obtained  from  the  plateau 

Steady-state  Acid  pump  OFF 


Fig.  4.  Transient  pressure  signal  of  pressure  transducer  P0  for  test  series  Al;  acid 
pump  is  switched  on  at  time  =  0  s. 


regions  is  judged  to  be  representative  of  the  reactor  during  sus¬ 
tained  operation.  In  the  following  discussion,  all  reported  values 
for  pressure,  flow  rate,  or  molar  concentration  are  with  reference 
to  the  average  value  recorded  over  the  steady-state  time  interval. 
The  molar  concentrations  of  the  acetic  acid  and  sodium  bicarbon¬ 
ate  (Cin>  add  and  Cm>  bicarb)  are  derived  from  the  time  averaged  flow 
rate  of  their  respective  pumps,  and  a  mass  balance  with  the  known 
solution  concentrations: 

Cin,acid  =  $ acid^ acid /V total  (2) 

C  in, bicarb  =  $ bicarb^ bicarb  /M bicarb^ total  (3) 

where  5add  is  the  acid  concentration  by  volume  (“strength”),  Sbicarb 
is  the  mass  loading  (g/L)  of  sodium  bicarbonate  and  Mhicarh  is  the 
molar  weight  of  sodium  bicarbonate. 

While  the  low  flow  rates  make  it  difficult  to  immediately  dis¬ 
cern  changes  in  the  axial  pressure  drop  due  to  the  reaction,  there 
is  a  prominent  rise  in  the  average  absolute  pressure  in  the  channel 
with  increasing  total  flow  rate.  There  are  two  possible  causes  for 
the  back  pressure  rise:  increased  hydraulic  losses  in  the  waste  lines 
and  valve,  or  expansion  of  gas  products  within  the  reactor.  Since 
the  hydraulic  diameter  of  the  tubing  following  the  reactor  outlet 
is  much  larger  than  that  of  the  reactor  channels,  increased  down¬ 
stream  pressure  losses  appear  less  likely  to  be  responsible  for  this 
phenomenon  than  the  sudden  evolution  of  gas  inside  the  reactor. 
This  is  also  corroborated  by  the  sudden  increase  in  reactor  pressure 
that  occurred  immediately  after  the  commencement  of  the  reac¬ 
tion.  While  the  molar  ratios  are  the  same  for  each  flow  rate,  greater 
amounts  of  gas  are  produced  by  the  neutralization  reaction  as  the 
molar  flow  rate  of  both  reacting  species  increases. 

This  is  further  confirmed  by  plotting  the  outlet  pressure  (from 
transducer  P4)  as  a  function  of  the  molar  flow  rate  of  the  limiting 
reactant  (i.e.,  the  reactant  with  the  lower  inlet  concentration  of  the 
two)  defined  as  V  ■  Cin>min  =  V  ■  mm(CinAdd,  Cin_bicarh),  where  V  is  the 
total  flow  rate  and  Cin>min  is  the  inlet  concentration  of  the  limiting 
reactant.  The  quantity  V  •  Cin  min  corresponds  to  the  maximum  flux 
of  C02,  in  mol/m2/s,  that  would  be  produced  assuming  that  the 
reaction  completed  within  the  channel.  Such  a  plot  including  all 
tests  is  demonstrated  in  Fig.  5. 


V'C  inmin  (mol/m2/s) 


Fig.  5.  Time-averaged  pressure  reading  at  pressure  transducer  P4  as  a  function  of 
maximum  molar  flux  of  C02. 
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Although  we  could  not  verify  the  remaining  quantity  of  unre¬ 
acted  species  at  the  reactor  outlet,  the  assumption  of  (approxi¬ 
mately)  full  conversion  within  the  reactor  was  checked  by 
measuring  the  reaction  times  in  a  stirred  beaker,  which  were  typ¬ 
ically  2  or  3  s;  the  residence  times  achieved  in  this  study  were  in 
the  range  of  1.5-3  s  based  on  a  reactor  volume  of  101  mL 
(±1  mL).  Additionally,  Sahoo  and  Sata  reported  reaction  half-times 
that  were  less  than  one  second  [7],  and  that  the  kinetics  were  lim¬ 
ited  by  mixing  time.  Because  of  the  high  degree  of  agitation  occur¬ 
ring  in  the  reactor,  and  the  residence  times  close  to  the  observed 
release  of  gasses  in  a  stirred  beaker,  it  can  be  concluded  that  the 
reaction  is  most  likely  very  close  to  completion  upon  reaching 
the  reactor  outlet. 

Comparing  the  change  in  the  values  of  back  pressure  rise  at  dif¬ 
ferent  flow  rates,  it  was  clear  that  the  frictional  and  gravitational 
pressure  drop  are  not  significant  factors  determining  the  reactor 
conditions.  However,  it  is  well  understood  that  the  prediction  of 
thermal  performance  of  heat  exchangers  is  closely  tied  to  the  pre¬ 
diction  of  frictional  pressure  drop  and  therefore  the  utility  of  the 
frictional  pressure  drop  measurements  in  predicting  reactor  ther¬ 
mal  performance  will  be  discussed  in  Part  II  of  this  study. 

The  absolute  pressure  data  obtained  from  the  outlet  was  also 
used  to  calculate  the  superficial  velocity  of  the  gas  phase  by  using 
the  ideal  gas  equation  of  state  to  obtain  the  density,  and  the  mea¬ 
sured  molar  flow  rate  (assuming  complete  conversion)  to  obtain 
the  mass  flux.  Fig.  6  illustrates  a  plot  of  superficial  gas  velocity 
(t4,m/s)  with  respect  to  gas-phase  mass  flux  (m",kg/m3).  While 
Ug  increases  with  greater  mass  fluxes,  diydm ”  appears  to  dimin¬ 
ish.  The  decreasing  slope  corresponds  to  an  increase  in  the  calcu¬ 
lated  gas  density  which  arises  from  the  enhancement  of  the 
reactor  pressure.  With  increase  in  mass  flux  values,  the  tortuous 
channel  geometry  and  low  hydraulic  diameter  increasingly  impede 
the  escape  of  the  gas  (reaction  products)  from  the  heat  exchanger. 
This  further  implies  that  the  gas  escape  kinetics  govern  the  reactor 
pressure. 

3.2.  Flow  visualization 

The  high-speed  camera  was  moved  to  three  fixed  axial  locations 
while  performing  the  experiments  in  a  given  test  series  (i.e.,  fixed 
molar  ratio  and  flow  rate).  This  enabled  the  monitoring  of  the  evo¬ 
lution  of  the  two-phase  flow  in  the  axial  direction  with  a  framing 
rate  of  1000  fps;  the  field  of  view  was  approximately 
13  x  26  mm.  The  camera  had  512  x  512  pixel  resolution,  which 
was  reduced  to  256  x  512  to  reduce  the  size  of  the  video  files. 
The  similarity  between  the  flow  patterns  observed  in  the  experi¬ 
ments  obviates  the  need  for  an  exhaustive  discussion  of  each  video 


B/W  PRINT 


Fig.  6.  Superficial  velocity  of  C02  at  axial  position  of  pressure  transducer  P4  as  a 
function  of  maximum  gas  flux  at  reactor  outlet.  The  uncertainty  for  the  calculation 
of  Ug  obtained  with  the  Kline-McKlintock  method  [8]  was  approximately  9%. 


that  was  acquired.  Instead,  a  representative  video  file  (obtained 
from  test  series  B3)  was  used  for  the  flow  pattern  analysis. 

Over  the  range  of  flow  rates  and  molar  ratios  tested,  two  dis¬ 
tinct  flow  regimes  that  were  observed:  a  froth-like  flow  (Fig.  7a), 
and  a  dispersed  bubbly  flow  (Fig.  7b).  The  froth  flow  was  observed 
where  the  reaction  was  particularly  strong,  such  as  near  the  inlet 
port  after  the  acetic  acid  and  sodium  bicarbonate  solutions  were 
mixed.  The  size  of  the  bubbles  decreased  as  greater  reactant  con¬ 
centrations  were  used,  indicating  that  the  bubble  size  is  inversely 
correlated  with  the  gas  evolution  rate.  Dispersed  bubbly  flow 
was  typical  of  lower  molar  flow  rates,  or  in  regions  where  the 
reaction  had  ceased  due  to  complete  consumption  of  one  of  the 
reactants.  In  such  cases,  there  was  a  wider  range  of  bubble  sizes, 
but  all  were  still  much  smaller  than  the  furrow  width  (10  mm) 
or  the  maximum  channel  spacing  (4  mm).  The  difference  between 
these  two  types  of  flows  shows  that  the  additional  agitation 
produced  by  the  chemical  reaction  acts  to  inhibit  the  formation 
of  larger  bubbles.  The  fine,  uniformly  distributed  bubbles  in  froth 
flow  also  explain  the  decrease  in  the  gas  superficial  velocity 
discussed  previously:  as  the  flow  becomes  more  homogenous  in 
nature,  the  relative  velocities  between  the  gas  and  liquid  phase 
decreases.  The  violent  churning  and  mixing  produced  by  the  reac¬ 
tion  ‘pulverizes’  the  flow  such  that  it  behaves  more  as  a  single  fluid. 
As  larger  quantities  of  gas  are  evolved  in  the  channel,  increased 
collisions  between  bubbles  and  the  tortuous  channel  passages  also 
contribute  to  the  homogenizing  effect.  In  no  instance  were  any 
large  gas-continuous  flow  structures  observed,  which  also 
supports  the  aforesaid  assertion. 

The  position  of  individual  bubbles  could  be  tracked  on  a  frame- 
by-frame  basis  to  reconstruct  their  flow  paths.  Due  to  the  similar¬ 
ities  between  the  video  footage  obtained,  the  flow  paths  shown 
here  are  representative  of  the  results  obtained  in  this  study. 
Fig.  8  shows  the  paths  of  several  bubbles  tracked  over  the  course 
of  15-20  frames  from  various  times.  Bubbles  of  different  size  and 
shapes  can  be  observed  following  nearly  identical  pathways;  thus, 
there  is  a  definite  steady  flow  structure  that  exists  even  as  the 
chemical  reaction  progresses  to  completion.  Fig.  7a  and  b  corre¬ 
spond  to  the  screenshots  shown  in  Fig.  8a  and  b,  respectively. 


Fig.  7.  Screenshots  of  steady-state  operating  for  test  series  B3,  (b)  between  P0  and 
PI;  (a)  between  P2  and  P3.  Total  flow  rate:  0.043  L/s;  Cmin  =  0.18  [mol/L].  The  solid 
red  line  represents  the  crest  of  a  furrow  and  the  dotted  red  line  represents  a  trough. 
Blue  lines  indicate  the  direction  of  crossing  flow,  with  the  dotted  line  representing 
the  lower  furrow  and  the  solid  line  the  upper  furrow.  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of 
this  article.) 
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The  pathlines  for  froth  flow  compared  to  those  of  bubbly  flow  ap¬ 
pear  very  similar,  despite  the  differences  in  bubble  sizes. 

It  is  also  worth  noting  that  the  same  flow  components  discussed 
in  the  literature  [10]  for  single-phase  flow  also  appear  to  be  pres¬ 
ent  in  the  two-phase  flows  studied  here.  Namely,  the  ‘S’-shaped  or 
abruptly  curving  pathlines  indicate  a  bubble  following  a  sinusoidal 
wavy  flow  pattern  (i.e.,  alternating  between  upper  and  lower  cor¬ 
rugations  in  the  main  flow  direction),  while  others  follow  the 
direction  of  a  single  corrugation.  The  fact  that  the  same  basic  flow 
components  characteristic  of  single-phase  flow  in  chevron  PHEs 
also  exist  for  two-phase  flow  is  consistent  with  the  observations 
reported  by  Tribbe  and  Miiller-Steinhagen  [11]  for  bubbly  flow.  Be¬ 
cause  only  the  bubbles  closest  to  the  window  could  be  observed 
and  tracked,  the  flow  fields  shown  in  Fig.  8  do  not  include  bubbles 
from  the  entire  channel  depth.  It  may  be  argued,  based  on  symme¬ 
try  between  the  upper  and  lower  plates,  that  the  patterns  observed 
in  the  upper  corrugations  are  similar  to  those  in  the  lower 
corrugations. 

3.3.  Two-Phase  Multiplier  Assessment 

The  similarities  observed  in  the  flow  patterns  (bubbles  follow¬ 
ing  paths  similar  to  those  of  single-phase  flow)  recorded  in  this 
study  to  those  reported  by  Tribbe  and  Miiller-Steinhagen  [11]  in¬ 
vite  a  comparison  to  the  two-phase  multiplier  correlation  they  pro¬ 
posed  from  their  pressure  drop  measurements  for  air-water 
vertical  upward  flow  in  a  chevron  PHE.  In  fact,  these  authors  also 
tested  the  same  model  of  PHE  used  in  this  study  (Alpha  Laval 
M3)  [11,12].  The  two-phase  multiplier  is  defined  as: 


where  the  subscript  F  indicates  the  frictional  pressure  drop  compo¬ 
nent,  and  TP  and  SP  represent  two-phase  flow  and  its  equivalent 
(i.e.,  same  mass  flux)  single-phase  flow,  respectively.  The  correla¬ 
tion  due  to  Tribbe  and  Miiller-Steinhagen  [12]  is  as  follows: 

<Dtp  =  1.423  -  0.0724  lnX  +  1.031  /X  (5) 

where,  X  is  the  Lockhart-Martinelli  parameter, 


Fig.  8.  Pathlines  traced  from  test  series  B3,  (a)  between  P0  and  PI;  (b)  between  P2 
and  P3.  Each  colored  line  corresponds  to  a  particular  bubble’s  path.  The  solid  black 
lines  represent  the  crest  of  a  furrow  and  the  dotted  black  lines  represent  a 
trough.(For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  this  article.) 


(6) 


The  authors  noted  their  correlation  to  be  valid  over  the  range 
0.06  <  X  <  10.  When  applied  locally  in  a  “stepwise  manner”,  Tribbe 
and  Miiller-Steinhagen  were  able  to  use  Eq.  (5)  to  reproduce  the 
total  experimental  pressure  drop  within  4%  using  [12].  We  likewise 
attempted  to  reproduce  the  total  pressure  drop  measured  during 
the  experiments  obtained  from  the  reactor  using  Eq.  (5).  In  two- 
phase  flow,  the  total  pressure  drop  is  typically  expressed  in  terms 
of  a  gravity,  acceleration,  and  friction  component  [9]: 


A  P  =  APg  +  A  PA  +  APF 

(7) 

where, 

rz 4 

APC  =  /  g[(l  -  a )p,  +  ccpg]dz 

JZq 

(8) 

A PA  =  r  p,Uj/(  1  -  a )dz 

JZq 

(9) 

A  PF  =  [\Uptff/2Dh)dz 

Jzq 

(10) 

Ui  is  the  equivalent  single-phase  superficial  velocity  (m/s)  (i.e.,  as  if 
liquid  were  flowing  alone  at  the  same  total  mass  flux),  p  is  the  den¬ 
sity  (kg/m3),  a  is  the  void  fraction, /is  the  single-phase  friction  fac¬ 
tor,  and  g  is  the  gravitational  acceleration  (m/s2).  The  subscripts  / 
and  g  signify  the  liquid  phase  and  the  gas  phase,  respectively.  These 
integrals  were  evaluated  between  the  position  of  pressure  trans¬ 
ducer  P0,  positioned  slightly  forward  from  the  inlet  port  (axial  dis¬ 
tance  z0),  and  pressure  transducer  P4  just  before  the  outlet  port 
(position  z4);  a  total  length  of  approximately  275  mm  as  measured 
parallel  to  the  channel  chevron  spine.  To  estimate  the  void  fraction, 
Kreissig  and  Miiller-Steinhagen  [13]  recommended  a  correlation 
due  to  Rouhani  [14],  which  was  also  used  here: 

a  =  3L  { [i  +  0.12(1  -  x)]  (x/pg  +  (1  -  x)/Pl) 

+l'18/«M,Pb5)  [ga(p,  ~ Pg)]°25}  1  (11) 

In  this  equation,  the  total  mass  flux  and  quality,  respectively,  are: 

Kta,  =  ^[P,  +  MgCin(\-r,)]  (12) 

x=(\-t1)(pl/MgCin  +  'l-tlp  (13) 

where,  r\  =  C/Cin,  C  is  the  mass  concentration  of  the  limiting  reactant, 
Cin  is  the  value  of  C  at  the  inlet  port,  and  Mg  is  the  molecular  weight 
of  the  gas  phase  (44  kg/kmol  for  C02).  In  order  to  simplify  the  eval¬ 
uation  of  the  integrals  in  Eqs.  (8)-(10),  it  was  assumed  that  the  gas 
mass  fraction  was  constant  throughout  the  channel  length.  This 
assumption  yields  the  following  approximations: 


APg  “  [(1  -  oc) Pi  +  a ;pg]  (z4  -  z0) 

(14) 

APa  -  PiUflCi  -  a)-'  -  1] 

(15) 

Since  the  frictional  pressure  drop  was  insignificant  compared  to  the 
absolute  pressure,  the  density  was  evaluated  at  the  outlet  pressure 
P4  and  room  temperature  (24  °C)  per  the  ideal  gas  equation  of  state. 

As  an  initial  assessment  of  their  suitability  for  gas-evolving 
reacting  flows  in  chevron  plate  heat  exchangers,  Eqs.  (5)  and  (11) 
are  used  in  a  1-D  plug  flow  reactor  model  to  estimate  the  pressure 
drop  in  the  HEX  reactor,  and  compared  to  the  experimental  results. 
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The  steady  species  transport  process  for  a  dilute  aqueous  solution 
of  acetic  acid  inside  the  PHE  channel  may  be  approximately  de¬ 
scribed  by  the  following  relation: 

|U-(k/l/,)»f  (16) 

where  /<  is  a  rate  coefficient  that  will  depend  on  fluid  mixing  ( since  the 
reaction  is  diffusion  limited),  and  will  therefore  vary  with  Ui  and  tur¬ 
bulence  generated  by  gas  bubbles.  For  the  sake  of  simplifying  the  cal¬ 
culations,  k/Ui  =  k *  is  taken  to  be  a  constant  (assuming  k  ~  Ui).  It  is  also 
dependent  on  temperature,  but  the  experimental  conditions  are,  to  a 
very  good  approximation,  isothermal.  As  mentioned  before,  the  reac¬ 
tion  between  aqueous  sodium  bicarbonate  and  acetic  acid  is  first  or¬ 
der  with  respect  to  acetic  acid  [6]  when  the  release  of  C02  is  the  rate 
limiting  step  [ 7  ].  This  transport  equation  may  be  integrated  to  obtain: 

r](z)  =  exp{— k*z}  (17) 

The  assumption  in  this  model  is  that  acetic  acid  is  the  rate  limiting 
reactant.  Otherwise,  to  account  for  the  exhaustion  of  a  rate  limiting 
reactant  at  axial  location  z*,  a  “cut-off’  will  have  to  be  imposed  such 
that  rj{z  ^  z*)  =  (Cin  -  Cmin)/Qn,  where  Cmin  is  the  molar  concentra¬ 
tion  of  the  rate  limiting  reactant  at  the  inlet.  The  axial  concentration 
profile  is  used  with  the  differential  forms  of  Eqs.  (7)-(10)  to  predict 
the  incremental  pressure  drops  over  I<  control  volumes  of  length  Az, 
which  are  summed  to  calculate  the  total  pressure  drop.  Setting 
I<  =  100  was  sufficient  for  convergence.  It  was  also  determined  that 
setting  k*  =  20  1/s  gave  good  agreement  with  the  experimental  re¬ 
sults  -  changing  k*  to  higher  values  did  not  improve  the  agreement, 
which  supports  the  assumption  that  the  majority  of  the  reaction 
had  completed  by  the  first  pressure  transducer  (P0). 

Calculating  the  appropriate  value  of  the  Lockhart-Martinelli 
parameter  (Eq.  (6))  required  an  appropriate  single-phase  flow  cor¬ 
relation.  Water/propylene  glycol  mixtures  of  20%  and  31.5%  pro¬ 
pylene  glycol  by  volume  were  used  to  obtain  a  correlation  that 
represents  the  single-phase  behavior  of  the  chevron  PHE  with 
the  acrylic  window  installed.  For  the  range  8  <  Re  ^  70,  the  data 
could  be  represented  by  the  relation /=  24/Re  to  within  an  average 
deviation  of  ±23%  and  an  expected  deviation  of  ±64%  at  a  95%  con¬ 
fidence  level. 

It  was  found  that  the  frictional  pressure  drops  were  severely 
underestimated  by  Eq.  (5).  A  correction  factor  of  3.0,  which  is  de¬ 
rived  from  the  average  of  the  ratio  between  experimental  <DTP  val¬ 
ues  and  those  predicted  by  the  uncorrected  Eq.  (5),  was  introduced 
to  approximately  compensate  for  the  simplifying  assumptions 
inherent  in  Eq.  (17).  After  applying  the  correction  factor  to  Eq. 
(5),  the  total  pressure  drop  was  reproduced  with  errors  ranging 
from  0.7%  to  49%,  with  an  average  deviation  of  26%  at  95%  confi¬ 
dence.  Some  of  this  disagreement  may  possibly  be  attributed  to 
the  assumption  of  a  constant  value  of  k *  for  all  of  the  experiments 
and  the  neglect  of  coupling  between  mixing,  chemical  kinetics,  and 
gas-liquid  interaction.  Representative  results  comparing  the  axial 
pressure  distributions  measured  experimentally  and  predicted 
using  the  method  described  above  (with  the  correction  applied 
to  Eq.  (5))  are  shown  in  Fig.  9.  The  deviations  from  the  expected 
linear  axial  pressure  profile  in  Fig.  9  can  be  attributed  to  the  pres¬ 
ence  of  sudden  changes  in  channel  geometry  at  the  distribution 
zones  where  pressure  transducers  P0  and  P4  are  located.  The 
experimental  two-phase  multiplier  values  and  the  empirically  ad¬ 
justed  Eq.  (5)  are  compared  in  Fig.  10,  where  the  predicted  values 
can  be  seen  to  fall  within  the  spread  of  the  experimental  values. 

Established  two-phase  flow  correlations  [1 1,14],  with  appropri¬ 
ate  empirical  adjustment,  provided  acceptable  engineering  accu¬ 
racy  for  pressure  drop  calculations.  However,  the  fact  that  the 
current  two-phase  multiplier  was  significantly  higher  than  that 
of  air-water  flow  (about  9-fold)  points  to  the  need  to  consider 
the  coupling  between  chemical  kinetics  and  hydrodynamics  if  high 


Fig.  9.  Axial  pressure  profiles  from  test  series  D1  (dashed  line)  compared  to  model 
predictions  (solid  lines).  Colors  differentiate  tests  with  slight  differences  in  flow  rate 
and  pressure  from  the  nominal  set  points  to  show  that  the  overall  pressure  drop  is 
repeatable.  Markers  correspond  to  axial  location  of  absolute  pressure  transducers. 
(For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 


Fig.  10.  Comparison  of  experimentally  calculated  values  (blue  diamonds)  and 
predictions  Eq.  (5)  with  multiplier  adjustment  of  3.0  applied,  (green  triangles)  of 
OTP.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 

fidelity  calculations  are  required.  The  additional  fluid  agitation 
produced  by  the  sudden  evolution  and  expansion  of  gas  within 
the  flow  is  likely  to  contribute  to  considerably  higher  frictional 
losses  than  if  the  liquid  and  gas  phases  were  simply  injected  at 
the  inlet  as  in  references  [11,12].  Although  the  multiplier  correc¬ 
tion  to  the  two-phase  multiplier  correlation  is  unable  to  fully  ac¬ 
count  for  this  effect,  it  is  able  to  provide  accuracy  sufficient  for 
estimating  pressure  drop  values  over  the  range  of  Lockhart-Marti¬ 
nelli  parameter  values  (X)  considered  (between  8  and  14)  for  the 
reactor  (and  reactants)  used  in  this  study. 

4.  Conclusion 

In  this  study,  hydrodynamic  behavior  in  a  chevron  plate  heat 
exchanger  with  multi-phase  flow  with  chemical  reaction  (involv¬ 
ing  gas  bubble  evolution)  is  investigated  experimentally.  The 
chemical  reaction  between  acetic  acid  and  sodium  bicarbonate 
provided  fundamental  insights  into  the  reactor  hydrodynamics, 
and  guided  the  selection  of  hydrodynamic  correlations  for  gas-gen¬ 
erating  reacting  flows.  High-speed  flow  visualization  showed  that 
the  vigorous  reaction  has  a  homogenizing  effect  on  vertical  upward 
gas-liquid  flow,  preventing  the  formation  of  any  large  heteroge¬ 
neous  flow  structures,  which  supports  the  use  of  correlations  based 
on  the  assumption  of  homogenous  flow.  Using  the  correlation  pro¬ 
posed  by  Tribbe  and  Muller- Steinhagen  [11,12]  for  the  two-phase 
multiplier,  with  an  empirical  correction  factor  of  3.0,  along  with 
the  void  fraction  correlation  developed  by  Rouhani  [14],  the  total 
pressure  drop  could  be  reproduced  with  a  maximum  error  of 
49%.  The  level  of  accuracy  obtained  using  these  two  correlations 
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was  typical  for  general  two-phase  flow  correlations,  suggesting 
they  are  appropriate  bases  for  developing  empirically  adjusted 
relations  for  estimating  pressure  drop  and  void  fraction  in  PHEs 
with  gas-generating  reacting  flows. 

Disclosure  Statement 

To  our  knowledge,  there  are  no  conflicts  of  interest  between  the 
authors  or  their  respective  institutions  with  the  International  Jour¬ 
nal  of  Heat  and  Mass  Transfer 

Appendix  A.  Test  matrix 

The  following  table  lists  the  experimental  parameters  for  each 
test  set  and  constituent  series.  The  flow  rates  are  average  readings 
during  the  steady-state  ‘plateau’  time  interval. 

List  of  experimental  parameters  for  each  test  set  and  constituent 
series.  The  flow  rates  are  average  readings  during  the  steady-state 
‘plateau’  time  interval. 


Test  series 

Al 

A3 

A4 

A5 

A6 

A7 

Test  set  A 

Acid  strength 

5% 

5% 

5% 

5% 

2.35% 

2.35% 

Bicarb  loading 

20 

20 

20 

20 

25 

25 

(g/L) 

Total  flow  rate 

0.064 

0.072 

0.074 

0.065 

0.041 

0.032 

(L/s) 

Cin,  acid 

0.14 

0.13 

0.13 

0.14 

0.14 

0.14 

(mol/L) 

Cin,  bicarb 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

(mol/L) 

Test  set  B 

B1 

B2 

B3 

B4 

Acid  strength 

5% 

5% 

5% 

5% 

Bicarb  loading 

20 

20 

20 

20 

(g/L) 

Total  flow  rate 

0.058 

0.049 

0.042 

0.07 

(L/s) 

Cin,  acid 

0.18 

0.19 

0.2 

0.18 

(mol/L) 

Cin,  bicarb 

0.19 

0.19 

0.18 

0.19 

(mol/L) 

Test  set  C 

Cl 

C2 

C3 

Acid  strength 

5% 

5% 

5% 

Bicarb  loading 

20 

20 

20 

(g/L) 


Appendix  A  ( continued ) 


Test  series 

Al 

A3 

A4 

A5  A6  A7 

Total  flow  rate 

0.068 

0.055 

0.031 

(L/s) 

Cin,  acid 

0.35 

0.34 

0.35 

(mol/L) 

Cin,  bicarb 

0.14 

0.14 

0.14 

(mol/L) 

Test  set  D 

D1 

D2 

D3 

D4 

Acid  strength 

8% 

8% 

7% 

7% 

Bicarb  loading 

30 

30 

31.4 

31.4 

(g/L) 

Total  flow  rate 

0.063 

0.053 

0.061 

0.037 

(L/s) 

Cin,  acid 

0.32 

0.33 

0.37 

0.4 

(mol/L) 

Cin,  bicarb 

0.28 

0.27 

0.26 

0.25 

(mol/L) 
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